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Abstract 

The flow about slender, pointed bodies can be characterized by different 
states with angle of attack. At moderate-to-high angles of attack (a ~ 40°), 
a steady, asymmetric vortex pattern develops along the body, leading to a net 
lateral force. At higher angles of attack (a ^ 60°), the aft-end of the body devel- 
ops an unsteady von Karman shedding. As the angle of attack approaches 90°, 
the entire body length exhibits a time-dependent vortex shedding pattern. The 
current work uses three-dimensional, thin-layer Navier-Stokes simulations to 
investigate the physical mechanisms of asymmetric vortex shedding at a = 40° 
and a = 60°. The development of an asymmetric vortex pattern via a convec- 
tive instability mechanism is investigated using tip bumps, surface roughness, 
and tip curvature. It’s found that surface roughness simulations can incite an 
asymmetric vortex state at a = 60° which is consistent with the application 
of a tip bump, and the experimentally observed flowfield. The unsteady von 
Karman vortex shedding on the aft portion of the body is also well resolved. 
The use of surface roughness did not incite a flow asymmetry at a — 40° , and 
it was necessary to simulate tip curvature at this angle of attack in order to 
generate an asymmetric vortex state. 


1 Introduction 

It is well known that the wake behind a cylindrical body aligned perpendicular 
to the oncoming flow can develop an alternating pattern of strong vortices. This 
alternate shedding of vortices from the body results in an instantaneously asymmetric 


flow field, while in the mean* the flowfield remains symmetric. The alternate shedding 
of vortices from a circular cylinder is commonly referred to as von Karman vortex 
shedding, or a von Karman vortex street. If a slender, pointed body of circular cross- 
section, such as a cone-cylinder or ogive cylinder, is placed at an angle of attack 
relative to the oncoming flow, it’s possible to observe a steady asymmetric flow state 
relative to the lateral symmetry plane of the body. This introduces the counter- 
intuitive notion that a nominally axisymmetric body can develop an asymmetric 
flow, not just instantaneously as in the flow about a circular cylinder in crossflow, 
but in a time-averaged sense*. 

The flow about a pointed, slender body of revolution at angle of attack can thus 
be broken down into several (experimentally-observed) 2 9 regimes (cf. Fig. 1.). The 
angles-of-attack where these flow regimes transition is dependent upon a number of 
factors, including body geometry, Mach number, Reynolds number, etc. At relatively 
low angles of attack, roughly cr < 20°, the flowfield remains symmetric with respect 
to the lateral symmetry plane of the body, even if the flow is perturbed slightly. 
The flow in this regime may remain attached, or boundary layer crossflow separation 
can occur if the '>r ,<T ^o of attack ’S sufficient, leading a pair of vortm^s which are 
symmetric relative to the lateral symmetry plane on the leeward side of the body. At 
the opposite angle of attack extreme, 70° < a < 90° (again approximately), the flow 
about the length of the body exhibits von Karman vortex shedding. In between these 
extremes, it’s possible for a nominally symmetric body of revolution to develop a non- 
zero mean lateral force. In the low-moderate angle-of-attack range, 20° < a < 40°, 
the variation of mean lateral force with roll angle forms a continuous periodic curve. 
In this regime, as the angle of attack is increased, the maximum absolute value of 
the lateral force also increases. As the angle of attack is increased into the moderate- 
high range, 40° < a < 70°, the variation of lateral force with roll angle can become 
nearly a square-wave, and the case of a symmetric mean flow becomes very difficult 
to achieve. This can be referred to as a “bi-stable” regime, or a bi-stable state; the 
lateral force becomes independent of roll angle except for a “switching” from positive 
to negative (and vice versa) as a critical roll angle is passed. 

Many of the previous experimental investigations noted that small changes near 
the tip of the model had a large effect on the observed vortex shedding pattern. De- 
gani and Schiff, 10 Degani, 11 and Degani and Tobak 7-9 proposed that imperfections 
near the tip of the model could incite a convective instability in the flow, leading to the 
development of steady asymmetric vortex patterns at high angles of attack. A con- 
vective instability is a sub-class of more general instability responses (cf. Ref. 12). This 
is shown schematically in Fig. 2 as the response to an impulsive perturbation. The 
behavior in Fig. 2a is stable, in that the response decays in all directions with time. 
Fig. 2b shows an example of an absolutely unstable flow. Here the response to the 
initial disturbance grows in both the upstream ( x/t < 0) and downstream ( xft > 0) 

‘mean in this work refers to time-averaging. 

tp i s also possible for a nominally symmetric body to exhibit a steady asymmetric flow state in 
crossflow due to an asymmetric boundary layer transition to turbulence(cf. Ref. 1). 
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unsteady, symmetric mean flow 
(von Karman shedding) 


Figure 1: Schematic of flow states about slender pointed bodies. The shaded lines 
represent vortices rotating in opposite directions. 

directions. A convectively unstable behavior shows growth only in a downstream 
traveling packet. Given this response, the flow eventually reverts back to the state 
before the perturbation was applied. 

Computational Fluid Dynamics (CFD) simulations have been used to demonstrate 
the possibility of a convective instability mechanism for these class of flows. These 




(a) Stable (b) Absolutely Unstable (c) Convectively Unstable 

Figure 2: Schematic of responses to an impulsive perturbation. 

investigations have primarily emphasized using the controlled environment of a simu- 
lation to perform numerical experiments that cannot be performed in a wind tunnel. 
By perturbing the geometry near the tip of tin uody, Dcgani et al. ,13 and Let y Cv. 
al. 14 were able to demonstrate that the flowfield in the low-moderate angle-of-attack 
range behaves as if it possesses a convective instability. In the work of Degani et al., 
an asymmetric bump perturbation was applied to the tip of an ogive cylinder body 
until an asymmetric flow state developed, then removed, and the flow reverted to the 
symmetric state that existed before the perturbation was applied. Levy et al. noted 
that in the low-moderate angle of attack range increasing the size of the perturbation 
led to an increase in the flow asymmetry, consistent with the idea of a convective 
instability. Similar investigations in the moderate-high angle-of-attack regime 7 have 
indicated that the onset of von Karman vortex shedding in this regime may be due to 
an absolute flow instability, although the numerical results are more difficult to inter- 
pret in this angle-of-attack regime. Some researchers have been able to simulate an 
asymmetric mean flow in the moderate-high regime without perturbing the flowfield 
in any manner, 15,16 however these simulations used conical flow, or other simplifying 
assumptions, which do not allow for the possibility of a convective instability. 

In Ref. 20, the author found that ogive cylinders placed at 40° angle of attack 
would exhibit an unsteady variation in normal force at Rep = 80,000 and 200,000. 
These variations were occurring at harmonics of the Strouhal frequency predicted 
by the Independence Principle, i.e. St a = St^D sin a, and hence were incorrectly 
related to a possible von Karman shedding mechanism. Subsequent investigations 18 
of crossflow separation in two and three-dimensions, and at several Reynolds numbers 
using varying grid densities, showed that these unsteady variations were caused by 
an instability in the free shear layer shed from the body. In a physical flow, these 
free shear layer instabilities would eventually transition and lead to tuibulence in the 
wake behind the body, as examined experimentally by Roshko 25 and later by Bloor. 26 
In Ref. 20, it was also noted that the amplification of round-off errors was sufficient to 
load to a von Karman vortex shedding along the length an ogive cylinder at rv = 60° 
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and Rcu — 200,000, and no geometric or other How perturbations was required to 
initiate this state. 

The current work used CFD simulations to investigate the physics of vortex shed- 
ding about pointed slender bodies; both the steady, asymmetric vortex patterns at 
moderate angles of attack, and the unsteady, von Karman vortex shedding that oc- 
curs at high angles of attack. Since a large database of experimental results exists for 
steady asymmetric vortex shedding for Reynolds number Reo = 10, 000 — 200, 000, 
the current work is performed in this same Reynolds number range. While it has been 
established that imperfections near the tip of a pointed body can incite a convective 
instability in the flowfield leading to an asymmetric vortex pattern, the physical mech- 
anisms by which this convective instability affects the changes on the flowfield are still 
not well understood. Without a detailed understanding of the physical mechanisms 
involved, designers seeking to eliminate or harness the lateral control power of asym- 
metric vortex shedding are left to test and develop concepts based solely on physical 
intuition. The current work seeks to provide insight into the physical mechanisms in- 
volved by performing controlled numerical experiments comparing different methods 
of perturbing the flow; tip bumps, surface roughness, and tip curvature. 


2 Approach 

2.1 Numerical Algorithm 

This work deals with numerical modeling of two- and three-dimensional separated 
flowfields. In order to accurately simulate the physics of these flows, it is necessary 
to use the Navier-Stokes equations as a physical model. The complete Navier-Stokes 
equations are often simplified using the “thin-layer” approximation 17 when performing 
high Reynolds number flow simulations. Neglecting the crossflow contributions to the 
viscous stress tensor still provides satisfactory results for the computed mean drag, 
and von Karman shedding frequency at the flow conditions of interest in this work 
(cf. Ref. 18), and all of the simulations in the current work utilize this simplification 
to reduce the computational costs. 

In order to accurately resolve the viscous stress terms, a fine grid spacing is re- 
quired in each coordinate direction that the viscous stresses cannot be neglected. 
This grid spacing makes the use of implicit methods desirable when using finite- 
difference algorithms because of the stiffness of the Navier-Stokes equations. The 
implicit schemes used in this work are all based on performing a time-linearization 
of the non-linear flux vectors, and obtaining a so-called “delta formulation” of the 
governing equations. 10 The backward Euler implicit time-integration scheme is used 
throughout. This is a first-order accurate scheme in time, however the allowable 
numerical time-step due to stability considerations for these flows is much smaller 
than required to resolve the physics of the unsteady von Karman shedding, and a 
first-order scheme has been found to provide satisfactory results (cf. Refs. 18,20). 
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In order to numerically integrate the governing equations, the spatial derivatives 
must be replaced by suitable finite-difference formulas. A straight-forward approach 
is to replace the inviscid derivatives with second-order, central-difference operators, 
and the viscous derivatives with centered, second-order mid-point operators. Unfor- 
tunately, this produces a scheme which is moderately unstable in three dimensions. 21 
This weak instability can be overcome however with the addition of numerical dissipa- 
tion of an order that doesn’t interfere with the accuracy of any physical mechanisms. 
The explicit dissipation is scaled by the local flow properties to reduce it’s influence 
within the boundary layer, where the natural viscous dissipation of the Navier-Stokes 
equations' is suitable. 

This combination of central-differenced spatial terms and higher-order numerical 
dissipation leads to what is referred to here as a block form of the Beam and Warming 
algorithm. 19,21 The block formulation is preferred over the diagonalized algorithm of 
Pulliam and Chausee 22 due to an inherent asymmetry of the diagonalization pro- 
cess that can lead to spurious flow asymmetries for high-angle-of-attack simulations 
(cf. Ref. 23). 

2.2 Computational Grids 

Computations were performed about a tangent-ogive cylinder configuration with 
a 3.5 diameter forebody attached to a 7.0 diameter cylindrical aft section. This is the 
same configuration examined computationally by Degani et al., 7 9,11 in previous work 
by the author, 20 and experimentally by Lamont. 5,6 In previous simulations of vortex 
shedding about a circular cylinder in crossflow and a tangent-ogive cylinder at a = 60° 
and a = 80°, it was found that the von Karman vortex shedding frequency was not 
well predicted at Re D = 80,000 or Re D = 200,000 (cf. Ref. 20). Grid refinement 
studies 18 determined a minimum crossflow resolution (A(f> = 1.0°, where <p is the cir- 
cumferential spacing) for the prediction of the von Karman shedding frequency, mean 
drag, and crossflow separation location at these flow conditions. This recommended 
crossflow resolution is utilized in this work. This leads to a grid with a crossflow plane 
having 249 circumferential and 58 radial points for a full-body (periodic) configura- 
tion, and 69 points in the axial direction. A typical axial grid plane is shown in Fig. 3 
(every third grid line is shown for clarity). The outer boundary of the computational 
domain is located fifteen body diameters from the surface. At the outer boundary 
a characteristic boundary condition is applied, while at the downstream exit plane 
a zeroth-order extrapolation of the flow cjuantities is performed which is consistent 
with the highly convective nature of the flowfield. 
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Figure 3: Lateral symmetry plane of the tangent-ogive cylinder computational grid. 

3 Numerical Results 

3.1 Tip Bump 

As was described in Sec. 1, the flow about a pointed slender body at moderate-to- 
high angles of attack can develop a steady asymmetric vortex pattern which produces 
a non-zero mean lateral force on the body. Previous CFD studies, 11, 13,14 were able 
to simulate an asymmetric vortex pattern on a slender body by placing a bump on 
the body very close to the tip of the nose. In Ref. 20, a similar numerical simulation 
was performed at o; = 60° and Rep = 200,000. This simulation utilized a similar 
crossflow resolution (A<p = 3.0°) as the previous studies by Degani et al. and Levy 
et al. This computation will be used as a baseline for the current work. 

The computed time-histories of the lateral and normal forces are plotted in Fig. 4. 
Initially, the geometry was unperturbed, and a von Karman vortex shedding pattern 
develops along the length of the body. This can clearly be seen in the variation 
of lateral force. The computed vortex shedding frequency and the circumferential 
location of the crossflow separation are not in good agreement with experimental 
investigations of shedding from a two-dimensional cylinder (of. Ref. 18). At t = 266 a 
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hump of height /w = 0.0 Id, where d is the local diameter was placed at xj D = 0.06 
and 90° from the windward symmetry plane (of. Appendix A). After the simulation 
had progressed through several cycles, the perturbation was removed at t 424. 




(a) Lateral force (b) Normal force 

Figure 4: Tip bump perturbation force coefficient histories (Moo = 0.20, a = CO , 
Re D = 200,000). 

When the bump is placed on the nose the flow develops a non-zero mean lateral 
force, while still maintaining vortex shedding about the aft-end of the body. The 
presence of the bump also causes the mean normal force to increase. This is consistent 
with experimental observations (cf. Ref. 4), which note that the greater the lateral 
asymmetry, the greater the normal force. The computed flowfield without the bump 
present contains two main vortices on the leeward side of the body that convect at a 
small angle to the longitudinal axis of the body and alternate in time. Figure 5 shows 
an instantaneous snapshot of the unsteady streaklines 28 after the tip perturbation 
has been applied. In the presence of the bump, the tip vortex on the perturbation 
side (R) quickly curves away from the body, w r hile the vortex on the opposite side (L) 
remains attached for a short distance, then curves away from the body. This steady, 
asymmetric pattern of vortices produces the non-zero mean lateral force. This steady 
vortex pattern is consistent with those observed experimentally 7 ’ 29 and previously 
simulated numerically by Degani 13 using a similar bump perturbation. Downstream 
of these two steady, asymmetric vortices the flowfield remains unsteady with the tip 
bump in place. 

When the bump is removed, the flowfield reverts to a symmetric time-averaged 
flow. This is consistent with the presence of a convective instability mechanism, as 
has been proposed by Degani and Tobak.'~ 9 

3.2 Surface Roughness 

The placement of a fixed bump on the tip of the ogive cylinder is a useful tool for 
examining the response of the flowfield to a particular perturbation, however it does 


8 



Figure 5: Instantaneous snapshot (t = 423) of streaklines with tip bump present 
(Moo = 0.20, a — 60°, Re D = 200, 000). 

have limitations. A criticism of the method is that a relatively large perturbation 
(compared for example to the local radius of curvature) is placed where the flowfield 
is known to be most sensitive. In fact, at low angles of attack it can be necessary 
to use large bump heights to produce the same level of lateral force as observed 
experimentally (cf. Ref. 14). 

The use of a large surface perturbation to “trigger” the convective instability 
mechanism may effect this change via a different physical mechanism than occurs 
experimentally. The large perturbation can cause the boundary layer to separate at 
the bump, and become entrained in the outer, inviscid flowfield. This could bias the 
flow towards a crossflow separation driven mechanism for asymmetric vortex shedding. 

It was investigated whether surface roughness can produce the same asymmetric 
flow states as observed in both previous numerical simulations with fixed tip bumps, 
and experimentally. The use of surface roughness more accurately mimics the phys- 
ical geometry encountered with an experimental model, and is easier to simulate 
numerically. The modeling of a bump near the tip can cause problems for both 
grid generation and numerical stability, while small surface perturbations are easy to 
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generate and do not appreciably skew the existing grid cells. 

The flow about a tangent-ogive cylinder at o = 60° was computed using a crossflow 
grid resolution of Ar/> = 1.0°. The lateral and normal force histories for the simulation 
are plotted in Fig. 6. The simulation initially had a "clean geometry, and again the 
flowfield develops a von Karman vortex shedding before the perturbation is applied. 
The variation of lateral force (and normal force) is larger than in the previous tip 
bump simulation, which is indicative of stronger vortices, and the vortex shedding 
frequency is lower. Both of these observations are consistent with the increase in 
crossflow grid resolution in the current surface roughness simulation compared to 
the prior tip bump studies. The vortex shedding frequency follow's very closely the 
curve St a = St 2D s\na suggested by the Independence Principle. At t = 155.5 
a random distribution of surface roughness with A mai = 1.0 x 10 -4 was applied 
(cf. Appendix A). After the surface roughness is applied, the flowfield develops a 
non-zero mean lateral force in the positive, or pilots-right, direction. The sign of the 
lateral force (configuration of the asymmetric vortex pattern) would essentially be 
chosen at random with this method. 




(a) Lateral force (b) Normal force 

Figure 6: Surface roughness force coefficient histories (¥«, = 0.20, a = 60°, Re D = 

200 , 000 ). 

Figure 7 shows an instantaneous snapshot of the unsteady streaklines computed 
with surface roughness applied. As in the bump simulation, one tip vortex (L) lifts 
rapidly from the body, while the lifting of the tip vortex on the opposite side (R) is 
delayed. An interface between the space-fixed oblique shedding and the unsteady von 
Karman vortex shedding on the aft-end of the body can be seen. 

Another useful tool for examining the asymmetric flowfield, and it s response 
to perturbations, is the variation of sectional lateral force with both axial distance 
along the body and non-dimensional time. This results in a three-dimensional carpet 
plot with sectional lateral force being the dependent variable. The sectional lateral 
force variation is plotted for both the tip bump simulation and the surface roughness 
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Figure 7: Instantaneous snapshot (t = 205) of streaklines with surface roughness 
present ( = 0.20, a = 60°, Rep — 200,000). 

simulation in Figs. 8 and 9 respectively. The lateral force maps are similar for the 
two simulations. In both cases th e un steady vortex shedding is clearly visible on the 
aft-end of the body (x/D > 6). The frequency of the vortex shedding is higher in the 
tip bump simulation, which is consistent with its lower crossflow grid resolution. On 
the tangent portion of the nose, both simulations show a steady lateral force which 
originates as a perturbation from the nose tip. In the tip bump simulation the lateral 
force is negative (towards the pilots left), while in the surface roughness simulation 
the sign is reversed. Both simulati on s also show a high amplitude, high frequency 
variation of lateral force near the nose-cylinder junction {x/D = 3.5). In the surface 
roughness simulation this variation is difficult to see because it’s superimposed on 
the steady lateral force on the nose region. It’s felt that these high frequency peaks 
are related to the onset of vortex shedding on the cylindrical portion, and the high 
frequency shear-layer instabilities discussed in Sec. 1. 

It isn’t possible to see in the resolution of Figs. 8 and 9, but very close to the 
tip point, the flowfield is highly unsteady in the tip bump simulation, while in the 
surface roughness simulation the flow is steady (and asymmetric). This unsteadiness 
is caused by the tip bump forcing a crossflow separation. The complete paper will 
include a detailed comparison of the flow in a crossflow plane very near the tip of the 
body for both the tip bump and surface roughness simulations. 
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Figure 9: Variation of sectional lateral F( 
0.20, a = 60°, Re D = 200,000). 


The effect of surface roughness on the flow about the tangent-ogive cylinder geom- 
etry was also investigated at 40° angle of attack. Initially, the same roughness height 
of A mni = 1.0 x 10~ 4 was applied, however this did not generate an appreciable 
lateral force. The roughness height was increased to A mn . r = 0.001 everywhere, yet 
still it was not possible to generate an asymmetric flow state. Finally, a configuration 
was created which maintained a roughness height of A max = 0.01 on the nose of the 
body ( x/D < 1.0) and A m(JI = 0.001 elsewhere, yet this still did not generate an 
appreciable lateral force on the body at this angle of attack. 

3.3 Tip Displacement 

It was attempted to determine whether it was possible to generate an asymmetric 
flow at a = 40° without resorting to the use of a tip bump. Another geometry modifi- 
cation which is consistent with experimental models is an asymmetry of the sharp tip. 
Moskovitz et al. 30 reported that the tip of their model appeared curved when viewed 
under a microscope. A configuration was generated where the tip point was offset 
0.001 D in the y-direction from the lateral svmmetrv plane, while the remainder of the 
geometry was unmodified. This flow about this geometry was simulated at cv = 40°. 
The lateral force immediately began to diverge from zero after the tip was moved, 
and eventually the flow reached a steady asymmetric state. A similar effect was found 
at a = 60°. The effect of the tip displacement on the the lateral force evolution can 
be traced to the perturbation response being convected down the length of the body. 
Figure 10 shows the variation of sectional lateral force along the body with time. 
The response convects downstream creating an asymmetric vortex strength, with the 
asymmetry in the positive direction for 3 < x/D < 6 and in the negative direction 
for 6 < x/D < 10. The final distribution of sectional lateral force is in relatively 
good agreement with the experimentally observed modes (cf. Refs. 4,6,14). From the 
off-surface streamline pattern (cf. Fig. 11) it was noted that the tip displacement in 
this configuration did not cause an asymmetric oblique vortex shedding pattern to 
develop, rather the two tip vortices have different strengths on opposite sides of the 
body. This is a prelude to the asymmetric vortex pattern that can develop, and is 
consistent with observations of a continuous “sine-wave distribution” of lateral force 
with roll angle. It was attempted to cause the vortices to “lift-off” the body by 
increasing the tip displacement. This did lead to a larger flow asymmetry, but the 
vortices remained nearly parallel to the body longitudinal axis. 


4 Future Work 

It is desirable to determine whether it’s possible to generate vortex “lift-off" at 
cv - 40° without resorting to the use of tip bumps. It will be attempted to incite 
this vortex state using solely the displacement of the tip point, however instead of 
simple lateral translations of the tip, vertical translations in combination will also be 
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Figure 10: Variation of lateral force with tip displacement present (M*. = 0.20, 
a = 40°, Re D = 200,000). 



Figure 11: Off-surface streamlines with tip displacement present (Moo = 0.20, a - 
40°, Re D = 200,000). 

attempted. These simulations will be included in the full paper along with a thorough 
comparison of the flow structures seen very close to the tip region in all three methods 
used to generate a flow asymmetry. 

A Geometry Perturbations 

A.l Bump Geometry Definition 

The bump geometry is defined as a circular perturbation following Fig. 12. h ma x 
is the maximum bump height, and is placed at an axial location X(,. The x-axis is 
assumed to coincide with the longitudinal axis of the body. The magnitude of the 


14 




(a) Front view ( y-z plane) 


(b) Top view ( x-y plane) 


Figure 12: Bump geometry. 


bump radius, |r^|, is also specified. A value of |r&| = 0.5|r(x)| has been found to give 
good results for tip bumps. The longitudinal distribution of bump contours is created 
by defining a bump height distribution in x, and building planes of circular bumps, 
i.e. 


hmaz 



Where 4 is the half-length of the bump perturbation. A value of 4 = 2|r(x(,)| has 
given good results with tip bumps. This leaves the origin of the bump circle to be 
solved from 


|r(:r)| + h(x) = \r a \ + \r b \ 

In the current work the bump was always placed 90° from the windward symmetry 
plane. 

In this manner the existing grid points (a;) in the crossflow plane will be modified 
only if they are located within the region of the bump that protrudes from the body, 
i.e. 

if (|cc - r a | < |r*ft|) then 

x ~r 0 + r b 

fi 

A. 2 Surface Roughness Simulation 

The effects of surface roughness were modeled by randomly applying small per- 
turbations to the geometry surface, 

v = r(x)(l + A) - A wax < A < A m „ T 


The maximum possible deviation between successive points is thus 2A mnT . 
turbation was only applied to half the surface points, 


The per- 


A = 


0 

±A mnx rand[0, 1] 


if rand[(), 1] < 0.5 
otherwise 


wdiere rand[0,l] is a pseudo-random number between 0 and 1 inclusive. The pertur- 
bation is also chosen to be positive or negative at random. The random numbers in 
the above equation (including the ± test) are considered to be distinct. 
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